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In recent years, further improvements in the efficiency of Cu2ZnSn(S,Se)4 photovoltaic devices 
have been hampered due to several materials issues, including cation disorder. Cu2SnS3 is a promis¬ 
ing new absorber material that has attracted significant interest in recent years. However, similar to 
CZTS, Cu2SnS3 displays cation disorder. In this work, we develop synthetic techniques to control 
the disorder in Cu2SnS3 thin films. By manipulating the disorder in this material, we observe crys¬ 
tal structure changes and detect improvements in the majority carrier (hole) transport. However, 
when the minority carrier (electron) transport was investigated using optical pump terahertz probe 
spectroscopy, minimal differences were observed between the ordered and disordered Cu2SnS3- By 
combining these results with first-principles and Monte Carlo theoretical calculations, we are able 
to conclude that even ostensibly “ordered” Cu2SnS3 displays minority carrier transport properties 
corresponding to the disordered structure. The presence of extended planar defects in all samples, 
observed in TEM imaging, suggests that disorder is present even when it is not detectable using 
traditional structural characterization methods. The results of this study highlight some of the 
challenges to the further improvement of Cu2SnS3-based photovoltaics, and have implications for 
other disordered multinary semiconductors such as CZTS. 


Introduction 

The development of efficient and scalable photovoltaics 
is an important challenge in today’s energy landscape. 
Although highly efficient thin film photovoltaics have 
been commercially available for some time, these tech¬ 
nologies may be limited in their scalability by the use 
of rare or toxic elements [1-3]. Cu 2 SnS 3 is a promising 
thin film photovoltaic absorber material which uses only 
abundant and nontoxic elements. Interest in this com¬ 
pound has rapidly increased in the past several years, 
and device efficiencies of 4% and 6% have been reported 
using Cu 2 SnS 3 -based absorber materials and alloys, re¬ 
spectively [4, 5]. 

Besides the reported device efficiencies, previous theo¬ 
retical and experimental works have revealed promising 
attributes of Cu 2 SnS 3 for photovoltaic absorber applica¬ 
tions. First principles calculations showed a wide range 
of phase stability, high optical absorption coefficient, and 
lack of Fermi level pinning in this material [6]. Experi¬ 
mental reports of the band gap of Cu 2 SnS 3 range from 
0.9-1.35 eV, well suited for a single-junction PV device 
[7, 8]. Furthermore, alloying with Si or Ge on the Sn site 
can serve to increase the band gap. Research has demon¬ 
strated that Cu 2 SnS 3 can be synthesized by several po¬ 
tentially scalable techniques such as sputtering and solu¬ 
tion processing [5, 9, 10]. 

A variety of crystal structures have been reported for 
Cu 2 SnS 3 thin films, including cubic, tetragonal, mono¬ 
clinic, and triclinic [5, 10, 11]. First principles calcula¬ 
tions for this compound suggest that the different crystal 
structures result from Cu/Sn disorder on the cation sites 
of the zinc blende-based lattice. In Ref. [12], the authors 
calculate the structures resulting from varying degrees of 
random disorder, ranging from cubic (fully disordered), 


to tetragonal (partially disordered), to monoclinic (fully 
ordered). A recent theoretical work has suggested that 
this disorder takes the the form of compositional inho¬ 
mogeneities caused by entropy-driven clustering (rather 
than fully random cation disorder), and could lead to 
potential fluctuations that negatively affect the carrier 
transport in Cu 2 SnS 3 [13]. It has been demonstrated 
that annealing Cu 2 SnS 3 thin films at higher tempera¬ 
tures promotes a transformation from the tetragonal to 
the monoclinic structure [14]. A few other studies have 
considered annealing of Cu 2 SnS 3 , including an investiga¬ 
tion of the behavior of extrinsic oxygen defects in this 
material [15]. However, there has been no work on the 
behavior of intrinsic defects during the annealing process, 
or the effects of annealing on the electronic properties of 
Cu 2 SnS 3 . 

The reported hole concentrations in Cu 2 SnS 3 are often 
high, ranging from 10 17 -10 19 cm -3 [16, 17]. This level of 
doping is generally considered too high for a photovoltaic 
absorber material, as it leads to a tunneling-enhanced in¬ 
crease in recombination within the absorber layer and at 
the heterojunction interface. In our prior work, we es¬ 
tablished that both Cu-poor and S-poor conditions are 
necessary for low carrier concentrations [18]. By control¬ 
ling both the Cu and S chemical potentials during film 
growth, we were able to tune the carrier concentration 
over 3 orders of magnitude and achieve films with p-type 
doping of less than 10 18 cm -3 . Even with this progress, 
lower carrier concentrations are still desirable for device 
integration, and have proven difficult to achieve through 
control of the Cu and S chemical potentials during film 
growth. When the Cu-poor and S-poor requirements are 
considered in conjuction with the Cu-Sn-S chemical po¬ 
tential phase space [6], it can be determined that the 


2 


lowest carrier concentrations should be found in Cu 2 SnS 3 
that is in equilibrium with SnS. 

In this work, we equilibrated our as-deposited 
Cu 2 SnS 3 with SnS by annealing the films under an 
SnS atmosphere. We found that this equilibration with 
SnS caused a transformation from the cubic/tetragonal 
structure to the monoclinic structure, along with a de¬ 
crease in the hole concentration. We used optical pump 
terahertz probe spectroscopy to investigate the minority 
carrier transport in our Cu 2 SnS 3 films. Both samples 
had short electron decay times of 0.1-10 ps, and similarly 
high degrees of charge localization, characteristic of a 
disordered structure. When the samples were imaged 
with transmission electron microscopy, both the cubic 
and monoclinic films had high densities of stacking faults 
and/or twins. We conclude that the existence of planar 
defects results in local disorder even in the monoclinic 
Cu 2 SnS 3 , and has a negative impact on the electronic 
transport in this material. The findings in this work 
have implications not only for Cu 2 SnS 3 , but also for 
Cu 2 ZnSn(S,Se )4 (CZTS) and other multinary disordered 
semiconductors. 

Methods 

Film synthesis and characterization 

The films described in this paper were deposited using 
combinatorial RF sputtering from 50 mm diameter CU 2 S 
and SnS 2 targets on heated 50x50 mm glass substrates, 
under conditions described in Ref. [18]. Further informa¬ 
tion about combinatorial synthesis approaches as applied 
to other absorber materials can be found in Ref.s [1! - 
22]. Characterization was performed on small (12.5 mm 
x 12.5 mm) sections of the original films, which were de¬ 
termined to be uniform with regards to composition and 
morphology (due to the shallow initial compositional gra¬ 
dient and size of film sections). X-ray fluorescence spec¬ 
tra were collected using a Fischer XDV-SDD instrument 
to obtain both the Cu/Sn ratios and the thickness of 
the films. The X-ray diffraction patterns were collected 
using a 0 — 20 geometry with Cu K-a radiation and a 
proportional 2D detector (Bruker D8 Discover with Gen¬ 
eral Area Detector Diffraction System software). Raman 
spectra were collected using a Renishaw inVia confocal 
Raman microscope configured with 532 nm laser exci¬ 
tation at 5% power, an 1800 mm -1 grating and CCD 
array detector. To compensate for the small spot size of 
the Raman measurement (~10 /im), six measurements 
were taken on each sample and averaged to provide a 
final spectrum. 

We used optical pump terahertz probe spectroscopy 
(OPTP), also known as time resolved terahertz spec¬ 
troscopy (TRTS), to measure the decay of pump induced 
conductivity in the ps range as well as the complex 
charge carrier mobility. The OPTP spectrometer has 
been described earlier [23]. The setup has been modified 
recently to record the pump induced change in THz 
reflection configuration. The transients were measured 


at the point of maximum change in THz amplitude, 
which represents the averaged THz conductivity de¬ 
cay, while the mobilities were analyzed by fitting the 
recorded THz reflection spectra to an optical model 
using the transfer matrix method. We assumed that the 
decay in conductivity is caused by a decaying carrier 
concentration where the maximum of the transient 
corresponds to all initially induced carriers. While the 
measured mobility is an average over all excited charge 
carriers, the conductivity is dominated by the electron 
(minority carrier) properties if the electron mobility can 
be assumed to be larger than the hole mobility [24]. 

Annealing experiments 

SnS powder was synthesized using Sn (Alfa Aesar, 
99.9%) and S (Alfa Aesar, 99.5+%) powders in a stoichio¬ 
metric ratio. The reactants were sealed under static vac¬ 
uum (<5 mTorr) in a quartz ampoule, heated to 400°C, 
and held for 10 hrs. Then, the reactancts were heated 
to 900°C and held for 18 hrs, before cooling to room 
temperature. This resulted in a phase pure SnS pow¬ 
der. The SnS 2 powder was synthesized by ball-milling 
the SnS powder with excess S for 2 hrs in a SPEX 8000D 
Mixer/Mill. The ball milled mixture was sealed in a 
quartz ampoule under static vacuum and held at 425° C 
for 18-24 hrs, resulting in a phase pure SnS 2 powder. 

Annealing was performed by placing the samples in a 
quartz tube along with ^1 g of SnS or SnS 2 powder in 
an alumina boat. A thermocouple was inserted into the 
tube such that the thermocouple was positioned directly 
over the sample. Flowing argon gas was supplied to the 
quartz tube, and exited via a bubbler filled with mineral 
oil. The samples were heated in a tube furnace with a 5 
hour total ramp time to the desired temperature. After 
the desired hold time, the furnace was turned off and 
allowed to cool to room temperature. 


Computational methods 

Density of states (DOS) and projected density of states 
were calculated using the VASP code [25] . Because of the 
band gap problem of the semi-local density functional 
theory (DFT) and large cells necessary to describe dis¬ 
order in Cu 2 SnS 3 we performed non-self consistent HSE 
calculations [26]. The wave functions for HSE were gen¬ 
erated with DFT -\-U [27] using Perde-Burke-Ernzerhof 
[28, 29] exchange correlation functional and on-site po¬ 
tential U = 7 eV applied to Cu(d) states. 

To compare degree of charge localization in the ordered 
(SG=Cc) and disordered Cu 2 SnS 3 we also calculated in¬ 
verse participation ratio (IPR). IPR measures the inverse 
of the fraction of atoms over which a given state is delo¬ 
calized and takes the form 


IPR = 


(Efc ?) 2 


(i) 


where the sums run over N atoms in the unit cell and ci 
are atom-projected density of states. For instance, for a 
state that is delocalized over all atoms in the unit cell, 



3 


O Ar 


10 


co 

> 


o 

CJ 


0.1 


30 min 


O Ar+SnS 2 
■ Ar+SnS 

■ 2 hr 

■ 8 hr 

- I l 

r O 


l 


Changing 


o 

anneal atm 


| ■ 



■ 


■ 

Increasing^*" 


■ 

anneal temp 




Increasing 



anneal time 




■ 

_i_ 


_l_ 


300 400 

Anneal temperature (°C) 


no 


500 


n= 

o 


Fig. 1 : As expected from the Cu-Sn-S chemical potential 
phase space, carrier concentration reductions are largest when 
the Cu2SnS3 films are annealed in equilibrium with SnS 
(square markers). Increasing anneal temperature shows some 
effect on the hole concentration; increasing anneal time causes 
more dramatic reductions. 


I PR = 1; for a state delocalized over half of the atoms 
in the unit cell, I PR = 2. 

For the disordered Cu 2 SnS 3 , the DOS and IPR 
were averaged over four independent atomic structures 
generated using Metropolis Monte Carlo method with a 
local motif-based model Hamiltonian described in Ref. 

[-3]- 

Results and discussion 

We investigated a variety of annealing conditions to 
determine their effects on the transport properties of 
Cu 2 SnS 3 . The annealing temperature, atmosphere, and 
duration were independently explored to determine the 
effects on the carrier concentration, and these results 
are presented in Fig. 1. As expected from the chemical 
potential phase space, the lowest hole concentrations 
were achieved under an SnS/Ar atmosphere (square 
data points in Fig. 1) [18]. Comparing the results for 
30 min, 2 hrs, and 8 hrs (red, blue, and green data 
points, respectively), we found that the kinetics of the 
annealing process were relatively slow. We determined 
that hole concentrations of ^5-10 17 cm -3 could be 
achieved by annealing under an SnS/argon atmosphere 
at 500° for 2 hours. These anneal parameters were 
used to investigate the structural and morphological 
changes that resulted from annealing, as well as both 
the majority (holes) and minority (electrons) carrier 
transport. Although lower carrier concentrations were 
achieved by annealing for 8 hrs (1-10 17 cm -3 ), the length 
of the anneal would preclude this sample from any kind 
of industrial scale-up, and as such, we have excluded it 
from further investigation. 


Structural changes & majority carrier transport 

We assessed changes in the crystal structure of the 
Cu 2 SnS 3 films using X-ray diffraction (XRD) and Ra¬ 
man spectroscopy. Overall, we saw a change from a 
cubic/tetragonal as-deposited structure, to a monoclinic 
post-annealed structure. This indicated a change from 
a disordered cation sublattice (cubic/tetragonal) to an 
ordered sublattice (monoclinic). 

Prior to annealing, the XRD pattern displayed only 
one peak at ^28.5°, shown in Fig. 2a. With only this 
peak, it is difficult to assign a definitive crystal structure 
to this XRD pattern: this peak could correspond to the 
cubic, tetragonal, or monoclinic structures of Cu 2 SnS 3 . 
Thus, we used Raman spectroscopy to further probe the 
structure of the preferentially oriented films. The pre¬ 
anneal Raman spectrum, shown in Fig. 2b, showed two 
major peaks at 299 and 351 cm -1 , corresponding to the 
cubic Cu 2 SnS 3 structure. Secondary peaks at 317 and 
338 cm -1 could be ascribed to the tetragonal crystal 
structure [8]. 

Post-anneal, the XRD pattern showed peaks corre¬ 
sponding to the monoclinic Cu 2 SnS 3 structure. The 
Raman spectrum had two major peaks at 290 and 351 
cm -1 , which can be assigned to the monoclinic struc¬ 
ture [30]. Two secondary peaks can be seen at 315 
and 371 cm -1 . In Ref. [30] these peaks are ascribed 
to a Cu 2 Sn 3 S 7 secondary phase, with the corresponding 
cation ratio (determined by energy dispersive X-ray spec¬ 
troscopy) of Cu/(Cu+Sn)=0.47. However, the cation ra¬ 
tio of our film, as determined by X-ray fluorescence, was 
Cu/(Cu+Sn)=0.63, suggesting that the amounts of sec¬ 
ondary phase present are minimal. 

Concurrent with the change in crystal structure 
upon annealing, we observed an improvement in the 
majority carrier transport. The carrier concentration 
decreased from 1.3-10 19 cm -3 for the as-deposited 
sample to 8.0-10 17 cm -3 after annealing. Additionally, 
the Hall mobility increased from 0.56 cm 2 /Vs to 8.2 
cm 2 /Vs. Two possible explanations for the increase 
in hole mobility are a reduction in grain boundary 
scattering due to grain growth during annealing, or a 
reduction in ionized defect density. Measuring the Hall 
mobility of the pre-annealed samples as a function of 
temperature shows the expected increase in mobility at 
low temperatures, suggesting that the hole transport 
is not grain boundary limited (see Fig. SI). Thus, we 
conclude that the reduction in ionized defect density is 
responsible for the increased hole mobility. 

Terahertz spectroscopy investigation of minority carrier 
transport 

In order to probe the minority carrier transport in the 
as-deposited and annealed Cu 2 SnS 3 , we performed opti¬ 
cal pump terahertz probe spectroscopy (OPTP). In con¬ 
trast to the improvement observed in the hole (majority 
carrier) transport, the changes in electron (minority car¬ 
rier) transport after annealing were less pronounced. 
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Fig. 3 : (a) The as-deposited sample shows two decay processes 
in the measured reflectivity: a 0.3 ps decay and a 7 ps decay; 
the annealed sample shows only the 7 ps decay process, (b) 
When the excitation power is lowered, the 0.3 ps decay can 
be detected in the annealed sample, suggesting that this state 
can be saturated for high carrier injection. 


Fig. 2 : (a) After annealing, the Cu2SnS3 films display XRD 
patterns corresponding to the monolinic structure. Prior to 
annealing, structural determination is difficult using only one 
peak, (b) Raman spectroscopy allows for elucidation of the 
pre-anneal structure as a mixture of cubic and tetragonal 
Cu2SnS3. The post-annealed spectrum shows peaks corre¬ 
sponding the monoclinic structure, as expected from the XRD 
pattern. 


In the as-deposited sample, the measured reflectivity as 
a function of time showed two decay processes: one with 
a decay time of 0.3 ps, and a second with a decay time of 
7 ps (see Fig. 3a). When the reflectivity of the annealed 
Cu 2 SnS 3 is compared to that of the as-deposited sample, 
it is clear that the 0.3 ps decay process is absent, but 
that the 7 ps decay process remains. However, when the 
excitation power is lowered to 2-10 17 cm -3 (as measured 
in terms of excited carriers at the surface of the sample), 
the 0.3 ps decay process can be detected even in the an¬ 
nealed sample, as shown in Fig. 3b. This suggests that 
the trap or defect states responsible for the 0.3 ps decay 
can be saturated for high carrier injection. At one sun 
conditions (lower than the excitation powers used here), 
this decay process may still be dominant. 

For both the as-deposited and annealed samples, the 
imaginary component of the terahertz conductivity is 
negative (Fig. S2). This negative value indicates charge 
carrier localization, which can be modeled using the 


Drude-Smith fit for conductivity [31]. Fitting our com¬ 
plex conductivity data results in a determination of the 
relaxation scattering time r and a constant ci, which 
represents the persistence of velocity after the first scat¬ 
tering event. The constant c\ can vary from -1 to 0, 
with ci=0 representing Drude conductivity and c\=-l 
representing complete carrier localization. For the as- 
deposited and annealed Cu 2 SnS 3 , the calculated value of 
ci is -0.91 (calculated 5-10 ps after excitation), indicat¬ 
ing strong localization in both samples. The relaxation 
scattering times were 33 fs and 69 fs, for the as-deposited 
and annealed samples, respectively. 

The minority carrier DC mobility of the samples can 
also be calculated from the Drude-Smith model (with 
cj= 0 for DC mobility), if the effective mass (m*) is 
known or approximated. Here, we approximate m* 
as 0.2m e , the value calculated for CZTS [32]. We 
calculate the minority carrier mobility to be 26 cm 2 /Vs 
for the as-deposited Cu 2 SnS 3 , and 55 cm 2 /Vs for the 
annealed sample. This factor of two increase in the 
electron mobility is directly related to the increase in 
the relaxation scattering time r. As with the increase in 
hole mobility, this is likely due to a reduction in ionized 
defect density after annealing. 

Theoretical explanation of terahertz spectroscopy results 

To interpret our finding that both the ordered and dis¬ 
ordered samples displayed high degrees of charge local- 
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Pre-Anneal 

Post-Anneal 

Structure 

Crystal structure (XRD) 

Cubic, oriented 

Monoclinic 

Crystal structure (Raman) 

Cubic / tetragonal 

Monoclinic 

Crystal structure (TED) 

Cubic, oriented 

Monoclinic 

Morphology 

Grain shape 

Columnar 

Equiaxed 

Grain size 

~50 nm 

200-500 nm 

Defects 

High density of {111} 
planar defects 

Some planar V other in¬ 
tragrain defects 

Majority carrier transport 

Hole concentration 

MO 19 cm -3 

8-10 17 cm -3 

Hole mobility 

0.56 cm 2 /Vs 

8.2 cm 2 /Vs 

Minority carrier transport 

Electron decay times 

0.3 ps, 7 ps 

7 ps 

Momentum relaxation time 

33 fs 

69 fs 

Electron mobility 

26 cm 2 /Vs 

55 cm 2 /Vs 

Electron localization constant 

-0.91 

-0.91 


TABLE I: Summary of changes in crystal stucture, film morphology, and majority and minority carrier transport caused by 
annealing Cu 2 SnS 3 thin films. Significant changes were evident in the structure and morphology, and a concurrent improvement 
was observed in the majority carrier transport. However, the minority carrier transport was mostly unaffected by the annealing. 


ization, we turn to the calculated inverse participation 
ratios (IPR) for ordered and disordered Cu 2 SnS 3 . The 
IPR values represent the degree of charge localization: 
an IPR value of 1 indicates that charge is totally delocal¬ 
ized; an IPR value of 2 indicates that charge is localized 
on \ of the atoms present in the structure. As shown 
in Fig. 4a, the ordered Cu 2 SnS 3 has IPR ratios varying 
between 1-2 at the band edges (similar to values for tra¬ 
ditional semiconductors such as Si or GaAs). In contrast, 
the disordered structure has high IPR values, especially 
at the conduction band edge, indicating a high degree 
of charge localization. Thus, when charge localization is 
used as a metric, both of our experimental samples could 
be considered “disordered”, as both display a high degree 
of charge localization (ci). 

At this point in the interpretation of our results, it is 
necessary to summarize some of the recent important 
theoretical work regarding the nature of disorder in 
Cu 2 SnS 3 [13]. The Cu 2 SnS 3 structure can be understood 
as an assembly of S coordination motifs, in which each 
S anion is tetrahedrally coordinated by four cations. In 
the ground state, Cu 2 SnS 3 is made up of only S-Cu 2 Sn 2 
and S-Cu 3 Sn motifs; the S-CU 4 , S-CuSn 3 , and S-Sn 4 
motifs are higher energy motifs and are not present. In 
the monoclinic (ordered) structure, the S-Cu 2 Sn 2 motifs 
are uniformly distributed throughout the supercell used 
in the computation. However, when the disordered 
structure is obtained from Monte Carlo simulations, the 
S-Cu 2 Sn 2 motifs form nanometer-scale clusters within 
the supercell. This clustering is due to the presence of 
attractive entropic forces in the disordered structure. 




Fig. 4: (a) Calcuation of the inverse participation ratio (IPR) 
for the disordered and ordered structures shows a high de¬ 
gree of charge localization (denoted by a high IPR value) for 
the disordered structure, suggesting that both of our samples 
exhibit disorder, (b) The disordered structure results in sig¬ 
nificant band tailing and a mid-gap state, both of which would 
significantly impact the performance of a Cu 2 SnS 3 -based pho¬ 
tovoltaic device. 
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Thus, when we speak of disorder in Cu 2 SnS 3 , we refer 
to this clustering effect, rather than a random cation 
distribution that would normally be implied. 

Microscopy investigations 

To understand how the monoclinic, ostensibly “or¬ 
dered” sample could still exhibit properties suggesting 
a disordered structure (such as high c\ values and ps car¬ 
rier lifetimes), we used transmission electron microscopy 
(TEM) to examine the microstructure of the two sam¬ 
ples. As shown in Fig. 5a, the as-deposited sample shows 
columnar grains (^50 nm across) that span the thickness 
of the film, and are strongly preferentially oriented in the 
<111> growth direction. All grains have an extremely 
high density of {111} planar defects (stacking faults and 
twins), which is confirmed by the visible streaks in the 
transmission electron diffraction (TED) pattern (see Fig. 
5c). The TED patterns also confirm the zinc blende (cu¬ 
bic) crystal structure, as was previously determined from 
the XRD pattern and Raman spectrum. The annealed 
sample shows a markedly different grain structure (Fig. 
5b). The grains are larger (200-500 nm) and equiaxed, 
with no detectable preferential orientation. In this film, 
some of the grains still have a high density of planar de¬ 
fects, while other grains do not exhibit visible defects. 
We note that this sample was annealed for 2 hrs, and the 
8 hr anneal resulted in a lower carrier concentration. It 
is possible that the 8 hr annealed sample has a higher 
fraction of grains without visible defects, and that even 
longer annealing times (although impractical for PV ab¬ 
sorber fabriction) would further reduce the density of pla¬ 
nar defects. The TED pattern for the annealed sample 
is complex, suggesting a monoclinic structure (Fig. 5d). 

The effect of stacking faults in Cu 2 SnS 3 can be under¬ 
stood by again considering this structure in terms of S 
coordination motifs. In the lowest energy state, only two 
motifs are present: Cu 2 Sn 2 -S and Cu 3 Sn-S. Each motif 
is coordinated by 12 other motifs, and a coordination 
number can be used to denote by how many motifs of 
the same type a certain motif is surrounded. The defect- 
free monoclinic structure requires that the coordination 
numbers are set as 2 and 7 for the S-Cu 2 Sn 2 and S-Cu 3 Sn 
motifs, respectively. When a stacking fault occurs, it al¬ 
ters the motif coordination numbers, causing local motif 
clustering and disorder. This is true even in the mono- 
clinic structure; from the TEM, we see that some grains 
are truly “ordered” in that no stacking faults are visi¬ 
ble, but other grains still have a high density of stacking 
faults, and this is ultimately what dominates the carrier 
transport in these films. 

If we now consider the calculated electronic density of 
states (DOS) for ordered and disordered (i.e., clustered) 
Cu 2 SnS 3 (Fig. 4b), the reason for poor minority carrier 
transport (Fig. 3) becomes clear. As a result of the high 
density of stacking faults (Fig. 5), both samples can 
be viewed as “disordered” when determining electronic 
properties. The DOS for disordered Cu 2 SnS 3 (Fig. 
4b) shows significant band tailing and a mid-gap state, 



Fig. 5: (a) As-deposited films show columnar grains with an 
extremely high density of {111} oriented stacking faults, (b) 
In the annealed films, the planar defect density is reduced in 
some grains, but remains high in others, (c) and (d) Tran- 
mission electron diffraction patterns confirm the zinc blende 
structure of the unannealed sample, and the complex TED 
pattern of the annealed sample suggests a monoclinic struc¬ 
ture. 


and both of these could contribute to the picosecond 
electron decay times observed in the OPTP (Fig. 3). 
Additionally, the band tailing reduces the band gap of 
Cu 2 SnS 3 (0.9-1.35 eV), ultimately reducing the final 
device performance. Moving forward, control of planar 
defects and local disorder will be an important research 
challenge that must be addressed for future improvement 
of Cu 2 SnS 3 photovoltaic devices. 

Implications of this work for CZTS and other disordered 
semiconductors 

Calcuations of the inverse participation ratio for 
Cu 2 ZnSnS 4 indicate similarly high degrees of charge lo¬ 
calization in the kesterite structure (see Fig. S3). The 
detrimental impacts of disorder in CZTS have been dis¬ 
cussed in multiple theoretical works, and include band 
tailing, potential fluctuations on the order of 200 meV, 
and the formation of nanoscale compositional inhomo¬ 
geneities [13, 33, 3^ ]. These challenges have proven dif¬ 
ficult to address experimentally, and may slow further 
increases in the efficiency of CZTS-based photovoltaics. 
Although Cu 2 SnS 3 is not nearly as well studied as CZTS, 
the similarities between these two materials suggest that 
Cu 2 SnS 3 -based devices may face similar development 
challenges in the future. 

Newer photovoltaic materials with potentially advan¬ 
tageous disorder effects are also being investigated, such 
as disorder-related band gap tuning in ZnSnN 2 [35, 36]. 




7 


These design ideas may offer advantages for the optical 
properties of the materials; however, the investigations 
will also need to take into account the effects of disorder 
on the carrier transport properties. Successful future 
investigations of disordered semiconductors will require 
the use of new modeling techniques to accurately assess 
the electronic effects of disorder, such as those developed 
in Ref. [13]. Experimental quantification of the effects 
of disorder is equally important, using not only struc¬ 
tural analysis, but also techniques such as the terahertz 
spectroscopy and high-resolution TEM used in this study. 

Conclusions 

In this study, we investigated the impacts of cation 
disorder on the electronic and structural properties of 
C^SnSs. We demonstrated the transformation from a 
cubic to a monoclinic crystal structure upon equilibra¬ 
tion of the Cu 2 SnS 3 with SnS, and a concurrent reduc¬ 
tion in hole concentration by almost 2 orders of magni¬ 
tude. We analyzed the as-deposited (cubic/disordered) 
and annealed (monoclinic/ordered) samples using optical 
pump terahertz probe spectroscopy, which detected 0.3 
and 7 ps decay processes in both samples. By analyzing 
the complex electrical conductivity from the OPTP mea¬ 
surement, we determined that both samples displayed a 
high degree of charge localization. When theory was used 
to calculate the charge localization in ordered and disor¬ 
dered Cu 2 SnS 3 , it suggested that charge localization is 
only found in the disordered structure. TEM investi¬ 
gations revealed high densities of stacking faults and/or 
twins in both samples, which could be responsible for 
local disorder and the associated poor minority carrier 
transport. 


Overall, the results presented in this work identify sev¬ 
eral challenges to the use of Cu 2 SnS 3 as a photovoltaic 
absorber. It is possible that different synthesis tech¬ 
niques may result in material without the high density of 
planar defects exhibited in our sputtered films. However, 
it is also likely that disorder in Cu 2 SnS 3 is unavoidable 
due to fundamental thermodynamic reasons (especially 
when the weak temperature dependence of the disorder 
is considered [13]). In this case, further efforts to 
improve upon the leading 4% efficient Cu 2 SnS 3 may be 
hampered by both the existance of local disorder, and 
the difficulty of detecting this disorder. As research into 
Cu 2 SnS 3 moves forward, it will be critical to assess the 
charge localization and disorder effects, and to correlate 
these properties with photovoltaic device performance. 
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Supplementary Information 

Investigation of annealing conditions 

We investigated a variety of annealing conditions to 
determine their effects on the transport properties of 
Cu 2 SnS 3 . Our experimental set up allowed us to easily 
vary the annealing temperature, atmosphere, and dura¬ 
tion. Each of these parameters was investigated indepen¬ 
dently to determine the effects on the carrier concentra¬ 
tion. 

From our previous work, we know that the lowest dop¬ 
ing levels in Cu 2 SnS 3 are achieved in Cu-poor and S- 
poor environments [ L8] . Consulting the Cu-Sn-S chem¬ 
ical potential phase space in Ref. [6] we see that these 
conditions are achieved when Cu 2 SnS 3 is in equilibrium 
with SnS. To verify this theoretical result, we annealed 
the Cu 2 SnS 3 films under three different atmospheres: Ar 
only, SnS 2 +Ar and SnS+Ar (note that SnS sublimes con- 
gruently). As expected, the SnS atmosphere resulted in 
the largest reduction in carrier concentration (see Fig. 
la). The SnS 2 atmosphere resulted in moderate car¬ 
rier concentration reductions, which is likely because the 
atmosphere was more S-rich than the SnS anneal. Fi¬ 
nally, the Ar only atmosphere produced very small re¬ 
ductions in carrier concentration. This can be attributed 
to the fact that this atmosphere is Sn-poor (analogous to 
Cu-rich) compared to the other two atmospheres investi¬ 
gated. Concurrent with the decreases in carrier concen¬ 
tration, we observed a slight increase in mobility for the 
sample annealed in the SnS atmosphere. This increase in 
mobility is likely connected to a decrease in defect-related 
scattering caused by the carrier concentration reductions. 

To investigate the kinetics of annealing in the SnS at¬ 
mosphere, we annealed the Cu 2 SnS 3 films at tempera¬ 
tures of 300, 400, and 500° C, while keeping the anneal du¬ 
ration constant at 30 minutes. Above 500°C, significant 
material loss and film delamination was observed. As 
shown in Fig. lb, the 400°C and 500°C annealing tem¬ 
peratures showed significant reductions in carrier concen¬ 
tration as compared to the 300°C sample. It is possible 
that temperatures higher than 500° C could result in fur¬ 
ther carrier concentration decreases, but as mentioned 
above, this was not possible due to material losses at 
higher temperatures. Thus, we also investigated longer 
anneal durations at 500° C. We observed that increasing 
the anneal duration from 30 minutes to 2 hrs resulted in 
a lower carrier concentration, and an 8 hr anneal time 
resulted in further carrier concentration decreases. This 
indicates that the kinetics of the defect reduction is fairly 
slow. Although it is possible that further decreases would 
be observed at anneal times longer than 8 hrs, longer an¬ 
neals were not investigated in this study. 



Fig. SI: The increase in Hall mobility at low temperatures 
suggests that the hole transport is not grain boundary limited 
in the pre-annealed samples. 
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Fig. S2: Carrier localization is suggested by the negative 
imaginary component of the frequency-dependent terahertz 
conductivity. 
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Fig. S3: The electronic effects of disorder are even more pro¬ 
nounced in the kesterite structure of CZTS than in Cu 2 SnS 3 . 

































